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The generalized Lyapunov exponents describe the growth of the second 
moments for a particular solution of the quasi- ID Schroedinger equation 
with initial conditions on the left end. Their possible application in the 
Anderson transition theory became recently a subject for controversy in the 
literature. The approach to the problem of the second moments advanced 
by Markos et al (cond-mat/0402068) is shown to be trivially incorrect. The 
difference of approaches by Kuzovkov et al (cond-mat/0212036, 0501446) 
and the present author (cond-mat/0504557, 0512708) is discussed. 

Recently Markos et al have published a comment [1] on the paper by Kuzovkov et al [2] 
where a growth of the second moments for a particular solution of the quasi-lD Schroedinger 
equation was related with the problem of Anderson localization. It was stated in [2] that 
the Anderson transition is of the first order and exists not only for space dimensions d > 2 
[3] but also in the 2D case. These statements look wild and the authors of [1] are right in 
not believing them. They are also right in statement that the growth of the second moment 
of wave function does not mean the growth of its typical value, which is governed by the 
average logarithm. However, Markos et al [1] go further and claim that analysis of the 
second moments presented by Kuzovkov et al [2, 3] is qualitatively incorrect and cannot 
provide any evidence for the metallic phase. These conclusions contradict the recent papers 
by the present author [4, 5] and are shown below to be trivially incorrect. The difference 
of approaches suggested in [2, 3] and [4, 5] is also discussed. 

Consider the 2D Anderson model described by the discrete Schroedinger equation 

and interprete it as a recurrence relation in the variable n, which we accept as a longitudinal 
coordinate. Initial conditions are assumed to be fixed on the left end of the system, while 
the periodic boundary conditions are accepted in the transverse direction, ipn,m+L = V'n.m- 
Cite energies Vn^m are considered as uncorrelated random quantities with the first two 
moments 

( Ki,m ) = , ( K(,mK',m' ) = W^^^n.n'^m.m' • (2) 

The growth of the second moments for this problem can be studied using the old idea by 
Thouless [6] based on the observation that variables ipn,m are statistically independent of 
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Vn,m with the same n. The main quantity of interest is {'ipnm)'i solving (1) for il)n+i,m and 
averaging its square, we can relate it with the pair correlators containing lower values of 
n. Deriving analogous equations for the pair correlators, we end with the closed system of 
difference equations for the quantities 

•^m,m' (^) = ( '4^n,m'^n,m' ) i 

ym,m'{r)) = (^n,m^n-l,m' ) , (3) 
^m,m' (^) = ( "071— l,m'0n,m' ) ; 

which for = has a form [4] 

Xm,m'{ri+1) = iy^5„i,„i/x^,TO/(n)+Xm+i,m'+i(n)+x„i_i,„i/+i(n)+x,„+i,^/_i(n)+x^_i,^'_i(n)+ 

+Xm,m'{n - 1) + ym+l,m'{n) + |/m-l,m'(^) + ^m,m'+l(^) + Zm,m'-l{n) (4) 
ym,m'{n+l) = -Xm+l,m'{n)-Xm-l,m'{n)-Zm,m'{n) 
Zm,m'{n+1) = -Xm,m'+l{n)-Xm,m'-l{n)-ym,ni'{n) . 

Instead to follow this natural procedure, Markos et al [1] invent their own approach. They 
rewrite (1) using the transfer matrix and construct the tensor product of two such matrices. 
After averaging, they arrive to a linear system of equations determined by a matrix 

f WH^l + Do^Do -Do 1 



10Do -1(2)1 
V 1(2)1 

where Dq — E — Hq and Hq is the Hamiltonian of the n-th slice for a pure system, 1 is the 
unit matrix of the size L x L. Even if being correct, this system is practically untractable 
due to sophisticated matrix constructions. The eigenvalues A = exp(ig) of the matrix (5) 
are declared to be determined by equation 

2 cos 2q - 2KiKj cos q + + - 2) = 2W'^i sin q (6) 

where Hi = E — 2 cospi and pi are allowed values of the transverse momentum p. Equation 
(6) has no resemblance with the corresponding Eq. 45 in [4]. The main difference is the 
absense of functions with the argument qL, which inevitably arise due to the boundary 
conditions and can be absent only in trivial cases. One can suspect, that this difference 
is related with incorrect treatment of the disorder term (like x 1 in (5)), without 

which the problem is trivial. This term is local inm — m' (see (4)) and nondiagonal in the 
momentum representation, which is seemingly used in (6). It looks likely that the local 
nature of this term was neglected and it was replaced by a suitable constant. We can try 
such thing for the system (4), replacing 5m,m' by unity. Then (4) is solved trivially 

xm,m'{n) = xe'P^^'P'^' e^" , ym,m'{n) = ye'^^^^^'^' e^" , Zm,m'{n) = ^e^f-+^^''-' e^" , 

(7) 



-l®Do 1(2) 

-101 







1 \ 



(5) 



2 



where allowed values for p and p' {2113/ L, s = 0,1,. — 1) are determined by the 
boundary conditions and the quantities x, y, z satisfy the equation 



/ + 4cospcosp' — 2sinh/? 2cosp 2 cos 7/ 
2 cosp 6^ 1 

2cos|j' 1 e^^ 




0. 



(8) 



The determinant vanishes under condition 



2 cosh 2/3 - 2kk' cosh /3 + (k^ + k'^ - 2) - 2W^ sinh /3 



(9) 



where K = — 2cosp, k' = — 2cosp'. In the case E = (when Eq.4 holds) this equation is 
identical to (6), if correspondence (3 = iq is taken into account. We see that, indeed, the 
disorder term in (5) was treated inadequately and its local nature was neglected. Physically, 
the equations (6, 9) correspond not to the true Anderson model but to its degenerate version 
when cite energies Kt,m are independent of m. 

In fact, the error is present already in the matrix (5). To obtain the term with disorder, 
one needs to produce averaging of the kind 
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and the result cannot be represented as W"^ 101. The latter form is valid for Vi = V2 in 
Eq.lO, while in the general case it corresponds to the model (1) with Vn^m being independent 
of m, in accordance with the previous analysis. Such model is of no interest and all 
conclusions made in [1] are irrelevant for the problem under consideration. 

We see that Markos et al [1] became victims of their own formahsm and were unable to 
make a way through the jungles of the tensor algebra. 



Now let us discuss the difference of approaches suggested by Kuzovkov et al [2, 3] and 
the present author [4, 5]. The initial system of equations (4) and its higher dimensional 
analogue (Eq. 5 in [5]) coincide with those used in [2, 3]. However, the quantity -2mm' (^) 
was not introduced in [2, 3] and its role was played by |/m'm(^)- As a result, the system of 
equations had no complete difference form and could not be solved in the natural manner 
with evaluation of full spectrum of exponents (3s- The Z-transform exploited in [2, 3] allow 
to find a solution only in the thermodynamical limit L — > cxd, using a questionable procedure 
of averaging over translations in the transversal direction. Translational invariance takes no 
place for the solution of (4) (see Eq. 10 in [5]) and the latter procedure probably eliminates 
all terms with the transverse momentum different from p = and p = G/2, if effectively 
only squares of these terms are relevant (G is a vector of a reciprocal lattice corresponding 
to the main diagonal of the Brillioun zone). As a result, a zero value for the critical 
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disorder was obtained for c? = 3 in the band center = [3], in a striking contrast with 
[5]. Fortunately for the authors of [2, 3], a condition p = G/2 corresponds to the minimal 
exponent (5min for d = 2 and d > 4, so the critical values ctq (corresponding to Wc in 
[4, 5]) were found correctly for these cases. As for the second critical point ctq for higher 
dimensionalities, we see no evidence for it in the spectrum of f3s- It looks that a filter 
function H{z) used in [2, 3] has not only poles corresponding to eigenvalues of the transfer 
matrix but also another singilarities, which are physically irrelevant. Correspondingly, we 
see no evidence of a special role of dimensionality d = 6, which is surely absent in the exact 
field theory approach [7]. 

Relation of the generalized exponents f3s with the Anderson transition was established 
in [4, 5] using the conventional variant of finite-size scaling [8]. Contrary, the papers [2, 3] 
were formulated in the engineer's language (using the concepts of signals, filters etc.), which 
has no direct relation to the Anderson transition. The limit L — > oo was taken in [2, 3] from 
the very beginning and the finite-size scaling approach could not be used for interpretation 
of results. As a consequence, a transition in the 2D case was interpreted as being of the 
first order, in evident contradiction with all available information. In fact, this transition 
is of the Kosterlitz-Thouless type [4] and there is no need for a serious revision in the 
weak localization region. The possibility of power-law localization is mentioned in both 
approaches [2, 3] and [4, 5], but in completely different contexts. 

Contrary to [2, 3], a clear distinction is made in [4, 5] between the generalized exponents 
j3s and the true Lyapunov exponents 7s : the latter are self- averaging quantities and surely 
have a more fundamental character. Fortunately, the knowledge of j3s provides essential 
information on ^s'- (a) inequality (3s > can be rigorously proven; (b) the order of 
magnitude relation ~ 7s takes place on the physical level of rigorousness; (c) /Sg and 7^ 
are practically equivalent from viewpoint of finite-size scaling philosophy. Relation of j3s 
and 7s with the parameters of the log-normal distribution is also of great importance. In 
fact, inequality (5s > '^Js is sufficient for the most responsible statements, such as existence 
of the 2D phase transition and absense of one-parameter scaling for 7^. Relation Ps ^ Is is 
used with great precaution and only in the cases when it does not contradict to numerical 
results. 

Finally, we do not consider as indisputable the conventional variant of finite-size scaling 
based on relation of the Anderson transition with the minimal Lyapunov exponent 7min- 
In the general case, some effective exponent 7e// should be used instead ^min [4]. Such 
modification restores one-parameter scaling in the weak localization region and eliminates 
the 2D transition from roughly half of models. 

This work is partially supported by RFBR (grant 03-02-17519). 
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